Coupling of Two Motor Proteins: a New Motor Can Move Faster 
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We study the effect of a coupling between two motor domains in highly-processive 
motor protein complexes. A simple stochastic discrete model, in which the two 
parts of the protein molecule interact through some energy potential, is presented. 
The exact analytical solutions for the dynamic properties of the combined motor 
species, such as the velocity and dispersion, are derived in terms of the properties 
of free individual motor domains and the interaction potential. It is shown that the 
coupling between the motor domains can create a more efficient motor protein that 
can move faster than individual particles. The results are applied to analyze the 
motion of helicase RecBCD molecules. 

PACS numbers: 82.39-k, 82.39Fk, 87.15Aa 

Motor proteins are active enzyme molecules that are responsible for generation of forces 
and molecular transport in biological systems 

ELQ 

. They move along polar molecular tracks 
such as cytoskeletal filaments and DNA molecules, and this motion is powered by energy 
released from a hydrolysis of adenosine-triposphate (ATP) molecules or related compounds. 
However, the mechanisms of the chemical energy transformation into the mechanical work 
are not fully understood 

Crystal structures of motor proteins reveal that they can be viewed as complex systems 
consisting of many functional domains Q, [J. It is assumed that this complexity of motor 
proteins appeared during the evolution as a way to perform simultaneously many biological 
functions and to increase efficiency. A good example is a RecBCD enzyme that belongs to 
a class of helicase motor proteins 0,0]. It processes DNA ends resulting from the double- 
strand breaks and other defects jjj . This motor protein unwinds the DNA molecule into two 
separate strands, and then it digests them by moving at the same time along the DNA until 
reaching a special signal sequence of bases, called the recombinational hot-spot (Chi site) 
RecBCD is a hetero-trimer made of three proteins: RecB (134 kDa), RecC 



flame 



;i29 kDa), and RecD (67 kDa) 



ABB. 



Two subunits, the individual RecB and RecD 
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domains, have helicase activities, consume ATP and act on 3'- or 5'-ended DNA strands, 
respects flflQ. Meanwhile, the third snhnnit (the dotnain ReeC) has no ATPase 
activity, and it functions as a clamp preventing the dissociation of the enzyme complex from 
the track jlfij . and it also recognizes the recombinational hot-spots jjj. Recent experiments 
provided data on DNA unwinding rates by RecBCD enzyme as well as by its active subunits 
RecB and RecD at the single- molecule level jf], 0, Surprisingly, the speed of the 

)rotein complex is significantly faster than the unwinding rates of individual motor proteins 
^]. In addition, RecBCD has a very high processivity, it can move more than 20000 bases 
without detachment, while the individual RecB and RecD subunits cannot travel more than 
!00 hases These ohservatrons raise the g enera! fondamenta, q uestion important for 

all biological systems: how the interaction between the different subunits is optimized to 
produce highly-efficient and multi-functional molecular biological nanomachines? The work 
presented here aims to address this specific issue by developing a simple stochastic model of 
the motion of motor proteins consisting of two interacting subunits. 

We assume that the motor protein complex consists of two particles, as shown in Fig. 
1. Each sub unit can move only along its own one-dimensional track that corresponds to 
the motion of RecB and RecD domains on the separate DNA strands. The positions of the 
particle A on the upper lattice is given by the integer /, while m specifies the position of the 
lower domain B. Due to the link between the motor subunits only the limited number of the 
motor complex configurations has to be considered. In the simplest description, we assume 
that 3 configurations are possible, i.e., < \l — m\ < 1 (see Fig. 1). More configurations and 
different interactions can be easily accounted in the model. Our approach is related to the 
theoretical model of a helicase motion proposed by Betterton and Julicher [18]. They view 
the DNA unwinding as a result of interaction between the whole helicase and DNA fork (a 
junction between double-stranded and single-stranded segments), while we discuss the effect 
of the molecular structure and the internal coupling of the subunits on the motion of motor 
proteins. 

The dynamics of the system can be described by a set of transition rates for domains A 
and B for discrete steps forward and backward along the corresponding tracks, as illustrated 
in Fig. 1. These rates depend not only on the type of the subunit, but also on the specific 
configuration of the cluster. However, the dependence of rates on DNA sequence is neglected 
in this approximation. For the configurations (I — 1,1) the upper subunit A can only move 
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forward with the rate u al , while the lower particle B can only hop backward with the rate 
w b2 . see Fig. la. Similarly, in the configurations (I + 1,1) [Fig. lc] the domain A can 
only go back with the rate w a2 , while the domain B can advance with the rate u b \. In the 
configurations (/, I) [Fig. lb] both particles can move forward and backward with the rates 
u a 2 (u b2 ) and w a \ (w b i) for the upper (lower) subunit. 

The interaction between the protein domains is specified by the parameter e > defined 
as energy difference between the states (I ±1,1) and (I, I), respectively (Fig. 1). We assume 
that the configuration (1,1) is energetically most favored, while the configurations (/ ± 1,/) 
have a higher energy. The possible reason for this might be the internal stress and/or the 
work needed to break the bond between the bases in DNA. It allows us to obtain the detailed 
balance thermodynamic relations for the transition rates: 

M = ^L expi+£/kB T), (1) 
^ = ^exp(-e/k B T), (2) 

W j2 Wj 

with j = a or b. The rates Uj and Wj are the forward and backward transition rates for the 
subunit A (j = a) and B (j = b) in the case of no inter-domain interaction, i.e., e — 0. 

We introduce P(l,m;t) as the probability to find the system in the configuration (l,m) 
at time t. It can be determined by solving of a set of independent Master equations, 

dP(l-l,l;t) 



dt 

dP(l + l,l;t) 



= u b2 P(l - 1, 1 - 1; t) + w al P(l, I; t) - (u al + w b2 )P(l - 1, /; t), (3) 
u a2 P(l, I; t) + w bl P(l + 1, / + 1; t) - (u bl + w a2 )P(l + l,l;t). (4) 



dt 

The corresponding equation for P(l, I; t) is just a linear combination of two equations pre- 
sented above, and therefore it is not considered. In addition, the probabilities satisfy the 
normalization condition, 

+oo 

{ p {l, h t)+P(l-l, l; t) + P(l + 1, /; *)) = 1, (all t). (5) 

l=— oo 

The solutions of Master equations can be found be summing over all integers — oo < I < +oo. 
Define new functions 

+oo +oo +oo 

p (t) = h t), Pi(t) =J2p(1- i, i; t), p[(t) =Y,P(i + i, J; t). (6) 

Z=— oo l=— oo l=— oo 
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Then, using the conservation of probability, the steady-state distribution can be easily de- 
rived, 

p = iM Pi = /3M A' = «M (7) 

where Q = 1 + a + (3, and a and (3 are given by 

M«2 + W b i U b2 + 10 a i , . 

a= ■ , 0= ■ . (8) 

Ubl + W a2 U a i + W b2 

From the knowledge of the probability densities the dynamic properties of the motor 
complex, such as the mean velocity V and the dispersion D, assuming the size of the step 



is one, can be calculated 
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The velocity is given by 



V = o {( Ual ~ W b2)Pl + ( u a2 + U b2 ~ W al ~ W bl )P + (u b \ ~ W a2 )P[} , (9) 



that yields the following result, 

V = ^(u a2 + u b2 -aw a2 - (3w b2 ). (10) 

Similarly for the dispersion we obtain 

n l/l r , , | n ] (V + w a2 )(u a2 -aV) {V + w b2 ){u b2 -f3V) \ 
D = o \ o Fa2 + u b2 + aw a2 + (3w b2 \ ■ ■ > . 

II I 2 U bl + W a2 U a i + W b2 J 

(11) 

If the motor domains are identical (A = B), then the transition rates are related as 

U a \ = U bl = Ui, U a2 = U b2 = U 2 , W a i = W bl = W X , W a2 = W b2 = U]_. (12) 

Then the expressions for the mean velocity and dispersion can be written in a much simpler 
form, 

y _ 2( Ul u 2 - w lW2 ) _ ^ 
ui+w 2 + 2(u 2 + u>i) ' 

u x u 2 + w x w 2 - V 2 
D = -. (14) 

Ml + W 2 + 2{U 2 + Wi) 

Now consider the case of symmetric motor domains without interaction (e = 0). If the 
particles are not connected in the cluster and allowed to move freely along its tracks, then 
their dynamic properties are 

V = u-w, D = (u + w)/2, (15) 



5 



with U\ = u 2 — u and w\ = w 2 = u. However, the average velocity and dispersion of the 
motor protein cluster (without interaction) are different. From Eqs. (JT5j) and (fT3j) we obtain 

n £ =0) = 2 l_ < D(e = 0) < 10 

V 3' 3 - D ~ 27 K J 

These results show that in the case of no inter-domain coupling the speed of the combined 
cluster is smaller than the rates of the free particles, as expected, although the fluctuations 
are also decrease. 

Much more interesting is the situation when the two motor domains interact with each 
other (e > 0). Using the detailed balance conditions (JTJ) and (J2J) the transition rates can be 
written in the following form 

Ujl = Ujj 1 - ' 1 , w n = wjj- ' 1 (17) 



u j2 = Ujj-^ 2 , w j2 = wrf-^ 2 (18) 

where 7 = exp(e/kBT), and j = a or b. The coefficients 9ji determine how the interac- 
tion energy is distributed between the forward and backward transitions Q, 18, Q|. They 
are closely related to the load-distribution factors that have been used successfully in the 
stochastic single-particle models of motor proteins Q]. It is reasonable to approximate the 
distribution factors < 8ji < 1 as equal to each other since they describe the energy of the 
similar processes in the motion of individual motor domains. However, the situation of the 
state-dependent energy-distribution factors can also be considered. 

For the motor protein complex with the symmetric domains (A = B) the effect of inter- 
actions can be analyzed by looking at the ratio of the cluster velocity to the velocity of the 
free non-interacting particles, 

r ^vr 5T7 (19) 

The dependence of the relative velocity on the interaction energy for different 9 is shown in 
Fig. 2. The most interesting observation is that for small values of the energy distribution 
factors (9 < 0.23) there is a range of interaction energies when the velocity of the motor 
protein's complex is faster than the velocities of the free particles. It contradicts the naive 
intuitive expectations, but it can be understood by considering again Eqs. 0, (|17|) and 
(JTHJ). Large interaction energies increase the transitions rates Uj\ and Wj 2 , and lower the 
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rates Uj2 and Wji for j = a or b as indicated in Eqs. (|17|) and (JTSJ). At the same time the 
probabilities of non-vertical configurations Pi and P[ are exponentially decreasing functions 
of the interaction energy. Then each term in Eq. (jOJ) has a maximum at some specific value 
of the interaction energy. Thus the dependence of the relative velocity on the interaction 
energy is a result of two opposing factors: the acceleration of some forward transition rates 
is balanced by the decrease of the probabilities of the configurations from which the motion 
is possible. In the limit of very large interactions the cluster will not move since it could 
only be found in the vertical configurations. 

The expression for the relative dispersion of the motor protein cluster with symmetric 
domains (A = B) is given by 



D 2^- e . . 27 1 - 6 
r D = 77- — 7T, — 9(u, w; 7) 



D 2 + 7 ^ ' ' (/ 2 + 7 



/ \ 2" 

2uw 47 ( u — w 



(20) 



(u + w) 2 (2 + 7) 2 \u + w / 

where it can be shown that 0.5 < g(u,w\ r y) < 1. The interaction energy also changes the 
dispersion in the similar way as the velocity, however, the effect is smaller. It can be seen 
from Fig. 3 where the ratio of the relative dispersion to the relative velocity is plotted for 
different transition rates. For all values of the parameters the relative dispersion is always 
smaller than the relative velocity. It is interesting to note that there are situations when 
ry > 1 and td < 1, i-e., the motor complex moves faster but fluctuates less than the free 
subunits, making it an extremely efficient motor protein. 

We can now apply our model for the analysis of the motion of RecBCD helicases. It was 
shown experimentally that at 37 °C and 5 mM ATP the mutant RecBCD* (the domain RecD 
is nonfunctional) and the mutant RecB*CD (the domain RecB is nonfunctional) unwind the 
DNA with rates 73 and 300 nucleotides/s, correspondingly [8|. Electron microscopy and 
biochemical assay data indicate that the free helicases RecB and RecD move with the same 
speeds as the corresponding mutant RecBCD* and RecB*CD enzymes j^. It can be argued 
that the backward rates Wji are small Q|, since the backward steps are rarely seen in the 
experiments on helicases. Then the average speed of DNA unwinding by the RecBCD 
complex can be approximated as 

y ~ frg + u ^ l ' d (21) 

7 + {uju b +u b /u a y 

where u a = 73 and u b = 300 nucleotides/s are the velocities of the free RecB and RecD 
proteins. Assuming 6 ~ 0, we obtain that V = 370 nucleotides/s for the interaction energy 



e ~ 6 kgT, in agreement with experimentally observed values of the RecBCD velocity 
The predicted energy of inter-domain interaction is very reasonable since it is larger than 2 
ksT needed to break the bond between the base pair in DNA [4, 3], but it is also smaller 
than a strong covalent chemical bond energy. Note also, that the maximal possible speed of 
the motor protein complex cannot be larger than the sum of the velocities of the individual 
domains. Thus RecBCD is working with almost maximal possible efficiency 

Our theoretical model has been greatly stimulated by the experimental observations on 
helicases. However, we believe that the effect of the interaction between the domains is very 
common for many biological systems, especially for motor proteins. Recent experiments 
on KIF3A/B kinesins, that are heterodimeric processive proteins involved in intraflagellar 
transport and Golgi trafficking, suggest that the interaction between the motor domains is 



important for maximizing the performance of these enzymes [2l|. By comparing the proper- 
ties of mutant homodimeric KIF3A/A and KIF3B/B proteins with the wild- type molecules, 
it was shown that the independent sequential hand-over-hand model cannot explain the 
heterodimer velocity data. The theoretical model similar to the one discussed above can 
be developed to describe the coordination between motor heads by assuming the interac- 
tion between them. We believe that taking into account the interaction between different 
domains is critical in construction of more realistic models of protein dynamics. 

The presented theoretical approach neglects several features that might be important 
for understanding the mechanisms of biological transport phenomena. Our description of 
the protein dynamics and the biochemical transitions is rather oversimplified. For example, 
the existence of intermediate conformations and states is not taken into account jl^]. We 
also assumed that the transition rates depend only on the configuration of the cluster and 
independent of the specific nature of DNA bases. At the same time the sequence dependence 



is rather weak for most helicases 



0. 



]. In addition, the effects of DNA elasticity and inter- 
domain flexibility have not been considered. 

In summary, the effect the inter-domain coupling in the motor protein complexes has 
been discussed by developing a simple stochastic discrete model. It was shown, using the 
explicit analytic formulas for the velocity and dispersion, that the interaction between differ- 
ent subunits in the enzyme complex might accelerate the speed of the cluster, as compared 
with the velocities of the free moving domains, without the significant increase in the fluc- 
tuations. This effect is due to the fact that the energy of interaction favors the compact 
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vertical configurations of the cluster, and it influences the forward and backward transitions 
differently. The asymmetry in energy-distribution factors results in a more efficient dynamics 
of the motor protein complex. Our theoretical method is used successfully to describe the 
dynamic properties of RecBCD helicases, and it is also relevant for other motor protein 
systems. Thus we present a possible explicit mechanism of the interaction between different 
domains in complex biological systems. It is interesting to note that the expressions for 
V and D [see Eqs. (II OJ) and (jllj) ] are identical to those obtained for a stochastic periodic 
model in which a single particle can move along two parallel tracks [22]. It means that the 
motion of two coupled particles can be mapped into the dynamics of the single particle on 
parallel-chain lattices. It indicates also that this approach can be generalized to include 
more motor protein's configurations and domains. Finally, the model can be also extended 
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to describe the effect of external forces on protein dynamics 
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Figure Captions: 

Fig. 1. Schematic view of a motor protein that consists of two domains. Transition rates 
u a i and w a i with % = 1 or 2 describe the motion of the domain A (small circles), while 
for the second particle B (large circles) the transitions rates are uu and wu- Only three 
configurations are allowed: (a) (I — 1, I) with the energy of interaction e > 0; (b) (/, /) with 
e = 0; and (c) (I + 1, 1) with e > 0. 

Fig. 2. The relative velocity for the complex motor protein with symmetric domains as a 
function the interaction energy for different energy-distribution factors 9. 

Fig. 3. The ratio of relative dispersion and relative velocity for the motor protein complex 
with symmetric domains as a function of the interaction energy for different forward and 
backward rates. 
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Figure 1: E. Stukalin, H. Phillips, A. Kolomeisky, Physical Review Letters. 
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Figure 2: E.B. Stukalin, H. Phillips, A.B. Kolomeisky, Physical Review Letters. 
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Figure 3: E.B. Stukalin, H. Phillips, A.B. Kolomeisky, Physical Review Letters. 



